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[~^ ' Abstract 
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At ultra-relativistic energies the minijet production in heavy-ion collisions be- 
comes sensitive to semi-hard parton rescatterings in the initial stages of the process. 
As a result global characteristics of the event, like the initial minijet density, become 
rather insensitive on the infrared cutoff that separates hard and soft interactions. 
This allows to define a nearly parameter-free saturation cutoff at which the initial 
^ , conditions may be computed. As an application we study the centrality dependence 

of the charged particle multiplicity, which is compared with present RHIC data and 
predicted at higher energies. 
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1 Introduction 

In heavy- ion collisions the partonic degrees of freedom of the two interacting nuclei become 
more and more important as the center of mass energy of the collision increases. At 
some point the main particle production mechanism in the initial stage becomes the 
liberation from the nuclear wave functions of a great number of partons, also called minijet 
plasma. At ultra-relativistic energies the partonic density of the two nuclei is so high that 
perturbative methods on one hand [0-0| and semi-classical non-perturbative methods on 
the other become applicable to the computation of the initial conditions of the minijet 
plasma. It's successive evolution will possibly lead to thermalization of the system and 
to the transition to the quark-gluon plasma phase, whose formation and characteristics 
depend crucially on such initial conditions. Though the latter are not directly accessible 
experimentally, they can be related to final state observables, like the charged particle 
multiplicity and transverse energy, allowing a test of the proposed theoretical models. 

We can divide in general the models in three classes: i) two-component models 0,0, 
in which particle production is assumed to be decomposable into the sum of a soft and a 
hard part according to some cutoff po; H) saturation models |1|,|^,^,0], which exploit the 
high parton densities involved in the process; iii) "others" , like the Dual Parton Model |^ 
and hydrodynamic models [0. To distinguish between them, it has been proposed in |^ 
to study the centrality dependence of the charged particle multiplicity, since this allows 
to disentangle to some degree the dynamical and the geometrical effects. For a review of 
the results of the above models on the charged multiplicity see [p!0| . 

At very high energies the target parton densities experienced by projectile partons 
are so high that the probability for them to have more than one semi-hard scattering 
may become non negligible already at the BNL Relativistic Heavy-Ion Collider (RHIC). 
At such regimes the usual perturbative computation p|, obtained by eikonalization of 
the minijet cross-section, may become inadequate. Indeed, it takes into account only 
disconnected two-parton interactions located at different points in transverse space but 
neglects the rescatterings. With the help of a few simplifying hypotheses semi-hard parton 
rescatterings have been included in the interaction mechanism in |TT|,p!2|, and lead to 



sizeable effects already at RHIC energies [|T^-|l5[|. Based on these results, in this paper 



we propose a new saturation mechanism for semi-hard minijet production and use it in 
a two-component model to compute charged particle multiplicities at RHIC and at the 
CERN Large Hadron Collider (LHC). 



2 Initial conditions and saturation 

When rescatterings are included in the interaction of two nuclei of atomic numbers A and 
B, the average number of A nucleus minijets at fixed impact parameter b is given by ||1 1|| : 



Ar7(6) = f d^rdxrA{x,b- r) [l - e-^/'i-'-«(--')rB(x',r)j ^ ^2.1) 

and the average minijet initial multiplicity is obtained by summing the analogous contri- 
bution from the B nucleus, N"^^ = N'^-' + N^-' . For simplicity we omit the flavour indices 



and consider only gluon-gluon interactions in our formulae, the inclusion of quarks being 
straightforward. In the numerical computations both the gluons and the quarks have 
been included. In Eq. (p.l|), Ta = TA{r)G{x) is the nuclear parton distribution function 
of the A nucleus, TA_{r) is its nuclear thickness function, normalized to A, evaluated at 
a transverse coordinate r relative to the center of the nucleus and G{x) is the parton 
distribution function of a proton at a given fractional momentum x. For simplicity we 
omit the flavour indices, an is the pQCD gluon-gluon cross-section at leading order in 
the high energy limit, 

anixx) = -nal^T I 1 — ] 6(xx's — Apl)6(l — x)9(l — x') , 

2 Po V xx'sj 

where we included all the kinematic limits and po is the cut-off that discriminates between 
soft and semi-hard interactions. We also included explicitly the /c-factor, k, to take into 
account higher order corrections. Both the cross section and the parton distributions 
depend on a scale Q = Po, which we take equal to the cutoff. In the numerical compu- 
tations we will set A; = 2 and use Woods-Saxon thickness function and GRV98LO parton 
distribution functions |1^ . 



Eq. ( p. ID may be interpreted as the integral of the average density of projectile partons 
(at a given x and r) times the probability of having at least one semi-hard scattering 
against the target. The exponent in Eq. (p.l|) may be interpreted as the opacity of the 
target nucleus, being proportional to the total transverse area occupied by its partons at 
the resolution scale po- Two interesting limiting cases may be studied. At high values of 
Po the target has a small opacity and is seen by the incoming partons as a rather dilute 
system. As a consequence N"^^ ~ 2 J (Prdxdx'T a{x, h — r)a h{xx')T b{x') , and we recover 
the usual perturbative result 0. On the other hand, at low values of po the target opacity 
increases: the target is becoming black to the projectile partons. As a consequence, the 
probability of scattering at least once becomes so high that nearly every projectile parton 
scatters and the minijet multiplicity reaches a limiting value instead of diverging as it 
happens in the Eikonal computation. 

In the regime where the target is almost black the semi-hard interactions are extracting 
from the projectile nucleus wave-function all its partons, and even if we use a lower cutoff 
no more partons are there to be extracted, for this reason the minijet multiplicity tends 
to saturate [|13|, see also Fig. |l]a. We call saturation cutoff the value of po at which this 



happens, and will denote it as psat- Of course the validity of this picture is limited to the 
kinematic regions where the saturation cutoff is in the perturbative range, psat ^ ^qcd- 
To give a quantitative definition of the saturation cutoff we start by considering a central 
collision of two equal nuclei. We define the upper bound for the minijet multiplicity as 

N'^mib = 0) = hm NT = f d'rdxTAix, r) . (2.2) 

Taking a very large fc-factor corresponds, indeed, to the limit in which the target becomes 
completely black and the semi-hard interactions are effective in extracting all the partons 
from the projectile nucleus. The limiting procedure is needed in order to keep track of the 
kinematic limits. As it is easy to see, N"^^ ~ ^um^ therefore we can define the saturation 
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Figure 1: a) The minijet multiplicity, N"''^ (solid line), and its limiting value, N™^ (dotted line), in a 
central Au-Au collision as a function of the cutoff po at RHIC and LHC energies. The dashed lines are 
the saturated minijet multiplicities, -/V™^, with a saturation parameter c = 0.7 (long dashes) and c = 0.8 
(short dashes). The intercept of the dashed lines with the dotted lines determines the saturation cutoff. 
b) The saturation cutoff as a function of the impact parameter at RHIC and LHC energies. The shaded 
area is the region where we estimate that the saturation criteria ceases to be valid. 



cutoff as tlie value of po sucli tliat tlie minijet multiplicity becomes a substantial fraction 
of its limiting value: 



N'^'{PO=Psat)=cNZiPO=Psat] 



(2.3) 



where the saturation parameter, c, is a positive number smaller than one. Notice that 
Psat = Psat{Vs,c) is a fuuctiou also of the energy of the collision. From our discussion 
it is obvious that c must be close to one to let psat he in the region where N"^^ is satu- 
rating. However, to stay in the perturbative regime we cannot choose it too close to one 
since Psat^O as c-^1. Finally, we define the saturated minijet multiplicity as the average 
multiplicity evaluated at the saturation cutoff: 



NTd 



NTd{V~s,c)=N-^^{po 



Psat) 



(2.4) 



In our approach this number represents also the multiplicity of partons produced in the 
early stage of the heavy ion collision. 

In Fig. |l]a we show the minijet multiplicity and its limiting value as a function of 
the cutoff Po at RHIC and LHC energies. The rapidity density at ?7 = is computed 
by integrating Eqs. ( p.l|) and ( p.2|) over a pseudo-rapidity interval |?7| < 1, where we 
approximated rj ^ log{x^/s/pQ), and by dividing the result by a factor two. The dashed 
lines represent the saturated initial conditions computed with c = 0.7 and c = 0.8. We 
can see that at a given energy N^l, which is obtained as the intercept of the solid and 
dashed lines, is nearly independent of the saturation parameter as long as the latter is close 
enough to one. Indeed, both at RHIC and LHC energy we obtain approximately a 3% 
increase in the saturated multiplicity going from c = 0.7 and c = 0.8. Therefore, whereas 
c is an arbitrary parameter its actual choice doesn't affect strongly the determination of 
the initial conditions. The dependence of psat on ^/s was studied in [Q where it is shown 
that for central collisions the saturation criterion is applicable from RHIC energies on. 

Unless we use nuclear thickness functions with sharp edges, like the hard-sphere dis- 
tributions, by applying blindly the saturation criteria to non central collisions we would 



obtain an impact parameter independent bound on the minijet multiplicity. Indeed we 
would have NJ2^{b) = j d?rdxT A{x.,h — r) = J d'^rdxrA{x,r). In this way, by requiring 
saturation as in Eq. ( |2.3| ) we would be asking the semi-hard interactions to extract all 
the partons from the projectile nucleus even in a very peripheral region, which is clearly 
unphysical. A simple way to implement the collision geometry in the saturation criterion 
is to cut by hand the thickness functions outside a given radius Re of the order of the 
nuclear radius. However, the minijet multiplicity as a function of the centrality of the col- 
lision turns out to depend too strongly on the choice of Re except at very high centrality 
or very high energies 



To find a less arbitrary way of implementing the collision geometry we look at the 
Glauber model computation of the average number of nucleons which participate in the 
collision: 



^^ part 



(h)=ld^rTA{h-r)VB{r) + A^B, (2.5) 



^B 



where Vsir) = 1 — [1 — app{s)TB{r)/B] and app is the inelastic pp cross-section, which we 



take from [jT^. At y/s = 130, 200, 6000 GeV we have app = 39, 42, 75 mbarn, respectively. 



Vb is the probability that a projectile nucleus at a given transverse coordinate r has 
at least one inelastic interaction with the target nucleons. Then, we may require the 
saturation only for the fraction of projectile partons that belong to a participating nucleon, 
and define an effective nuclear distribution function 

Tab{x, b, r) = Ta{x, b - r)VB{r) ■ 

Correspondingly, we have an effective minijet multiplicity, A^^ (b) = f d'^rdxTAB{x, b, r) [l- 

exp (^ — J dx'aHixx')Tb{x'))~\ , and an effective upper limit, N AUmi^) = / d^rdxVAB{x, b, r), 
which is no more 6- independent. Then, the saturation criterion generalized to arbitrary 
impact parameter becomes: 

N'^^PO = Vsat) = cNuLiPO = Psat) ■ (2.6) 



Finally, having determined psat in this way we use it in Eq. ( p.4| ) to compute the average 
initial parton multiplicity. 

In Fig. |l|b we show the saturation cutoff as a function of the impact parameter at 
RHIC and LHC energies for different saturation parameters c. The horizontal line show 
the limit of approximately 0.7 GeV whose intersection with Psatib) sets the limit of validity 
of the present approach, as will be discussed in the next section. 

Notice that the saturation cutoff, and consequently the initial conditions, are practi- 
cally determined by the choice of the parton distribution functions. As explained above 
and in [^] the initial conditions are nearly independent of the remaining free parameters, 
namely the saturation parameter c and the fc-factor. 

In the proposed mechanism saturation is reached when there are no more partons 
that semi-hard interactions can extract from the nuclear wave functions. In this sense 
this mechanism is a saturation of the minijet production and is intermediate between 
initial and final state saturation. In initial state saturation [jl| (see also 0,0) the parton 
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Figure 2: The minijet average occupation number in the transverse area (sohd hues): n™"* = A™-'/^^", 
where A™-' — N™^ {po) x tt/pq and A^" is the transverse area of a gold nucleus. The dashed lines show the 
average occupation number of the saturated minijets as a function of the saturation momentum. When 
n^-' > 1 the minijets begin to overlap transversely. 



density inside the incoming nuclei saturate due to a compensation between parton splitting 
and parton fusion processes in the DGLAP evolution, which induces a corresponding 
saturation in the minijet multiplicity. On the contrary, in the final state mechanism ||^,^ 
the saturation is assumed to be caused by the high density of produced minijets, which 
screens softer parton production due to parton fusion processes in the final state. In 
particular these final state interactions are assumed to set in when the transverse area 
occupied by the minijets becomes comparable to the nuclear overlap area. Both processes 
may, therefore, complement our saturation mechanism since the former modify the input 
parton distribution functions and the latter deals with a later stage process. However, as 
we can see in Fig. |^, at a/s = 130 GeV and y/s = 200 GeV the saturated minijets fill 
the transverse area only partially. Therefore final state saturation effects should not alter 
significantly our computations at RHIC energy, but may play some role at LHC 



3 Charged particle multiplicity 



We want to apply the saturation criterion for the semi-hard parton production in the initial 
stage of the collision to the computation of the charged particle multiplicity. Thanks to 
the self-shadowing property of the semi-hard interactions |J15|,|18|], even if in Eq. (|2.1|) 
only the semi-hard cross-section, an, appears, we are actually taking into account all the 
partons that had at least one semi-hard scattering, while their other scatterings may be 
semi-hard or soft with no restrictions. Therefore we are missing only the purely soft part 
of the production mechanism. This leads us to adopt a two-component model in which 
the charged particle multiplicity is written as the sum of a soft and a semi-hard part: 
dN''^/dr]{h) = dNf^^^J d7]{h) + dN^\/dr]{b). The soft part is assumed to scale with the 
number of participants, Eq. (^.51), so that M 



dr] 



(b) 



Npart{b) 



(3.1) 



Here npp{s) the pseudo-rapidity density of charged particles produced at r] = in pp colli- 
sion at a given cm. energy y^. We use the fit |l^, npp(s) = 2.5— 0.25 log(s) -(-0.023 log^(s). 
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The coefficient x = x(s) is a parameter that allows to adjust the relative weight of soft 
and semi-hard interactions and will be determined from the experimental data. Further, 
we assume the semi-hard part to be completely computable from the saturation criterion 
for minijet production described in the last section. To convert the minijet multiplicity 
to charged particle multiplicity, we further assume isentropic expansion of the initially 
produced minijet plasma and parton-hadron duality, so that 

where the factor 0.9 is due to the different number of degrees of freedom of the system 
in the minijet-plasma phase and in the hadronic phase ||^. To mark out the contribution 
of the hard part it is customary to divide the charged multiplicity by the number of 
participant pairs, so that the observables we are interested in are: 

-ih) = xn,,{s) + ____iA(5) (3.3) 



N,,rtib)/2 dv '^ '"" N,,rtib)/2 d^ 



^Kl. idN-'- 



and the fraction of semi-hard interactions, Fg.h. = i'^' j- ^ ■ 

To make a comparison with experimental data we have first to relate the observables 



appearing in Eq. ( |3.3| ), which are functions of the impact parameter, to the experimental 
ones, which are obtained as averages over centrality classes of events [20-2^. Following 0, 



23| , to which we refer for the details, we do this by studying the minimum bias multiplicity 
distribution of charged particles and by dividing the events in suitable subsets over which 
the average is performed. The next step is to extract the parameter x in Eq. (|3.1|) by 
comparing the computation for the 3% most central events and the phobos data at 
V^ = 130 GeV from Ref. [^. This value is then used to make predictions at higher 
energy. 

In Fig. 1^ we show both the results for the semi-hard part before the averaging over 
the centrality classes, and the results obtained after the averaging and the inclusion of the 
soft part. For each curve the result obtained by setting c = 0.7 and c = 0.8 in Eq. (|2.6|) 
is shown. 

At y/s = 130 GeV we find x = 0.445 and x = 0.453, for a saturation parameter c = 0.7 
and c = 0.8, respectively. These values of x correspond to a fraction of semi-hard interac- 
tions Fsh^ = 0.805 and F^./i, = 0.817, respectively, and show a good stability with respect 
to c. The relatively large value of F^./i. with respect to the common expectation of nearly 
a half and to the value of 0.37 extracted from phobos data in Ref. is due to the fact 
that we considered as belonging to the non-soft part of the observable also a semi-hard 
region 0.7 GeV ^ Po ^ 2 GeV. Note that we can push our perturbative computations to 
such low values of the cutoff because inclusion of parton rescatterings results in a rather 
small sensitivity of global observables on po in that region [|13], ^ . 

The two curves start with a moderate slope at high centrality and at some point they de- 
crease very fast. This happens when the corresponding saturation cutoff becomes smaller 
than 0.7 GeV, approximately. The reason for this behaviour is that the distribution func- 
tions are fitted just down to a scale Q ~ 0.9 GeV and they are numerically extrapolated 
at lower scales. Below 0.7 GeV the extrapolation gives an unnaturally fast decrease of 



cv 



1^ 18 

^ 16 

^ 14 

S 12 

10 



.] I 1 I I I I I I I I I I I I 1 I I 

_x piionos [20] i PIiriNIX [22] 
m PHOHOS [21] 



I I I I I I I I I I I I 



Vs = 300 GeV _ 




100 200 300 400 

— - c = .7 

c=.7 averaged 



c = .t 

c=.8 averaged 



100 200 300 400 

N 



part 



Figure 3: Charged particle multiplicity per participant pair, Eq. (3^), as a function of the number of 
participants at different center of mass energies and saturation parameter c = 0.7 and c = 0.8. In each 
panel the lower pair of curves (dashed and dotted lines) represent the semi-hard contribution. The upper 
pair of curves (solid and dot-dashed lines) are obtained by averaging the semi-hard contribution and by 
adding the soft part. At y/s = 130 GeV the parameter x is extracted from the 3% most central phobos 
events at ^/s = 130 GeV pi| . At higher energies the upper curves may be considered an upper bound, 
while the lower ones give a lower bound, see text. 



the parton densities, which results in the rapid fall of the minijet production. Then, we 
define the region of validity of our computations as the one such that Psat ^0.7 GeV, or 
in other words the one to the right of the knee in the charged multiplicity. 
The value of psat at fixed centrality decreases when the saturation parameter c increases 
(see Fig. |l|), therefore the curve with c = 0.8 is reliable for a smaller range of centrality 
than the curve with c = 0.7. They agree, however, in the common region of validity 
(showing a slight tendency to increase their slope when increasing c), and after the ex- 
perimental averaging and the fit to the most central data point, both describe well the 
experimental data. 

At -^i = 200 GeV we don't have any data to normalize the multiplicities to. However, 
the fraction of semi-hard to soft interactions is expected to grow with the energy of the 
collision, and we can use the value of Fg.h. determined at ^/s = 130 GeV to obtain an 
approximate upper bound for the charged multiplicities: for y/s > 130 



xn 



pp 



. . 1 -^■^■fe-lb=0,v/i=130 GeV C^A^g./i./c?^ /, 



«-'i-lb=0,^/i=130 GeV 



Npart/'^ 



0,sl 



(3.4) 



The curves for the two values of c agree over a wider range of neutralities. This is to 
be expected since the saturation cutoff at fixed centrality grows with the center of mass 
energy, and goes below the critical value of 0.7 GeV at smaller centrality. Notice also that 
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the slope of the curves has increased. 

At LHC energy, y/s = 6 TeV, the particle production is generally believed to be almost 
completely semi-hard. Therefore we expect that the data will be close to the averaged 
semi-hard multiplicity without any normalization (which is very similar to the lower curve 
plotted in Fig. 0). Though the saturation criterion is applicable over all the centrality 
range considered (see Fig. Po) the slope of the curves is rather sensitive to the saturation 
parameter, resulting in a larger theoretical uncertainty. We expect that a better treatment 
of the scale Q and of the pseudo-rapidity, which are taken to depend simply on the cutoff 
Po, could solve at least partially this problem. However, the average slope has increased 
confirming the trend observed at lower energies. 



4 Conclusions 

The inclusion of semi-hard parton rescatterings in the interaction dynamics of heavy- 
ion collisions at very high energy allows a reliable computation of the initial conditions, 
like the minijet multiplicity, and the introduction of a nearly parameter-free saturation 
criterion to determine the infrared cutoff to be used in the perturbative computations. 
The proposed saturation mechanism is intermediate between the initial and final state 
ones in that it deals with the saturation in the production of minijets. 

We tested our approach against RHIC data on the centrality dependence of charged 
multiplicities by using a two-component model in which the semi-hard part is assumed 
to be completely given by the proposed saturation criterion. At ^/s = 130 GeV we find a 
good agreement with the data, which allows us to extrapolate the results at the highest 
RHIC energy of ^/s = 200 GeV and at LHC energy, ^/s = 6 TeV, by putting upper 
and lower bounds on the charged multiplicities per participant pairs as a function of the 
number of participants and by predicting their slope. 
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